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Abstract 

We present a derivation of the stress field for an interacting quantum system 
within the framework of local density functional theory. The formulation is 
geometric in nature and exploits the relationship between the strain tensor 
field and Riemannian metric tensor field. Within this formulation, we demon- 
strate that the stress field is unique up to a single ambiguous parameter. The 
ambiguity is due to the non-unique dependence of the kinetic energy on the 
metric tensor. To illustrate this formalism, we compute the pressure field for 
two phases of solid molecular hydrogen. Furthermore, we demonstrate that 
qualitative results obtained by interpreting the hydrogen pressure field are 

not influenced by the presence of the kinetic ambiguity. 
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I. INTRODUCTION 



The stress, or the energetic response to deformation or strain, plays an important role 
in linking the physical properties of a material (e.g. strength, toughness) with the behavior 
of its microstructure. In addition, the spatial distribution of stress is an invaluable tool 
for continuum modeling of the response of materials. The stress concept has been applied 
at atomistic scales as well. Over the last fifteen years, there has been a continuing trend 
toward understanding various structural and quantum-mechanical phenomena in materials 
in terms of their response to stressB. 

For example, the residual stress at equilibrium has been used to assess the structural 
stability of systems containing surfaces or strained interfaces. It has been demonstrated 
that the desire to minimize surface stress can give rise to reconstructions on high symmetry 
surfaces! - !, and the stability of epitaxially grown bimetallic systems has been attributed to 
the formation of incommensurate overlayers, defects, and dislocations which minimize the 
stress near the metal-metal interfaced The stress can have significant effects on chemical 
reactivity as well. It has been shown that small molecule chemisorption energies and reaction 
barriers on certain strained metal and semiconductor surfaces are quite different from those 
on the corresponding unstrained surface00. 

Formally, studies of the above phenomena must include a quantum-mechanical descrip- 
tion of the system's electronic degrees of freedom. Therefore, one must consider how a stress 
is defined quantum mechanically. Methods for calculating the stress in quantum mechanical 
systems have been developed since the birth of quantum theory^!. However, research in de- 
veloping a formalism for determining the quantum stress in solid-state systems has recently 
been revitalized. This is mainly due to ever-increasing opportunities to perform accurate 
and efficient quantum-mechanical calculations on systems which exhibit stress mediated 
phenomena. 

The stress is a rank-two tensor quantity, usually taken to be symmetric and therefore 
torque-free. Two useful representations of the stress tensor are the volume-averaged or total 
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stress, Tij, and the spatially varying stress field cr^(x). The two representations are related, 
since the total stress for a particular region in a system is the stress field integrated over the 
volume. Nielsen and Martin developed a formalism for calculating the total quantum stress 
in periodic systems^. They define the total stress as the variation of the total ground-state 
energy with respect to a uniform scaling of the entire system. This uniform scaling corre- 
sponds to a homogeneous or averaged strain over the entire system. They further demon- 
strate that the total quantum stress is a unique and well-defined physical quantity. Their 
formulation has been successfully implemented to study a variety of solid state systems^. 
Other formalisms for determining the total quantum stress have been created as well& 

Although these formalisms have provided important tools for studying quantum stress, 
the stress field is a more useful quantity that contains important information regarding the 
distribution of the stress throughout the system. A knowledge of the spatial dependence of 
the quantum stress is vital if one wishes to predict the spatial extent of structural modifica- 
tions or understand phenomena at interfaces in complex heterogeneous systems. However, 
the quantum stress field can only be specified up to a so-called gauge term. Therefore the 
quantum stress field is not unique. A traditional way to develop a quantum stress field 
formalism is to consider the stress field's relationship with the force field. From this per- 
spective, the stress field can be defined as any rank-two tensor field whose divergence is the 
force field of the system: 

n*) = v^(x). (i) 

(Note that the Einstein summation convention for repeated indices is used throughout the 
paper.) One can add to o % i a term of the form 

where A l ^ k is any tensor field antisymmetric in j and k, and recover the same force field, 
thereby demonstrating the non-uniqueness of this stress field definition. General formu- 
lations for computing stress fields in quantum many-body systems have been derived by 



Nielsen and Martin, Folland, Ziesche and co-workers, and GodfreyEMlhlx Regardless of 
issues concerning non-uniqueness, the quantum stress field has been used to investigate ma- 
terial properties. For example, Ramer and co-workers developed a method to calculate the 
resultant stress field from an induced homogeneous strain in order to study piezoelectric 
effects in perovskites0. They incorporate the additional constraint that the field must be 
the smoothest fit to the ionic forces. However this method cannot be used to calculate the 
residual stress field at equilibrium, nor can it determine the energy dependence on strains 
which do not have the periodicity of the unit cell. Filippetti and Fiorentini developed a for- 
mulation of the stress field based on the energy density formalism of Chetty and Martini3ii. 
Since this formulation is based explicitly on the energy density, which cannot be uniquely 
specified, their stress field is also ambiguous. 

In order to ascertain the exact nature of these ambiguities and to determine their effect 
on using the quantum stress field as an interpretive tool, we have developed a formulation of 
the quantum stress field based on differential geometry. This approach is similar to Mistura's 
work in deriving pressure fields for classical inhomogeneous fluidsii. We first develop the 
necessary geometric framework and discuss the origin of ambiguities in the stress. We then 
derive the stress field in the context of local density functional theory (LDA-DFT). Finally, 
we apply our formulation by using the pressure field to explain the energetic ordering between 
two phases of solid molecular hydrogen. 

II. METRIC AND STRAIN FIELDS 

The geometric formulation for the stress field is developed from the relationship between 
a strain tensor field and the Riemannian metric tensor field. This relationship has been 
derived for flat manifolds^ and for manifolds of arbitrary curvature!!. The infinitesimal 
squared distance between two points on a Riemannian manifold can be written as ds 2 = 
gijdx l dx\ where the metric is a rank two symmetric tensor field. The system can then be 
infinitesimally distorted so that the positions of the points change, dx' % = dx % + du\ which 
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changes the squared distance to ds' 2 = gijdx H dx' 3 . However, a coordinate transformation can 
be made so that the new distance can be computed from the point of view of the predistorted 
positions: ds' 2 = g' i jdx i dx 3 , with the new metric defined as 

9m = 9ij (4 + dku*) (Sf + «V) • (3) 

The strain field e^- can be related to the difference between distances before and after a 
distortion has been applied: ds' 2 — ds 2 = 2eijdx l dxK Therefore the metric g'^ can be related 
to the strain: 

-<U ■ .'/</• ( 4 ) 
We define the stress field a 13 via the virtual work relationship between energy and strain 

5E= f y/gj'SeijdPx, (5) 



where g is the determinant of g^. Using Eq. |, the stress field can be expressed as 

1 5E 2 5E 



y/9 feij y/9 &9: 



where g^ is the inverse of the metric. It is this definition which we will use to compute the 
stress at the quantum length scale, with E = (ty\H\ty) defined as the energy of the ground- 
state |^) of a system described by a Hamiltonian operator H. This method of deriving the 
stress field is identical to that used to compute the "improved" energy-momentum tensor in 
field theories which are coupled to a background gravitational source@. Specifically, we shall 
be interested in the value of the stress field evaluated at the Euclidean metric = 5ij). 

III. UNIQUENESS 



For the virtual work theorem, (and therefore Eq. ||), to specify a unique stress field, 
two criteria must be obeyed. First, the metric must be varied freely, without constraints, 
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when computing the functional derivative. Second the quantum-mechanical total energy 
must be unambiguously defined for all metrics. Before we can explain the origin and the 
consequences of these criteria, we first must discuss in more detail the relationship between 
a metric tensor and the geometry of the underlying space containing the physical system. 
Metric tensor fields, for our purposes, can be classified into two categories: Those fields which 
cause the Riemann curvature tensor to vanish at every point in space (e.g. the Euclidean 
metric = Sij), and those fields which give rise to a non-zero Riemann curvature tensor 
(e.g. the metric for the two-dimensional surface of a sphere). The latter class of metrics will 
be referred to as non-Euclidean. The Riemann curvature is a rank four tensor field defined 
as 

Rikim = 7^ {dkdigim + did m gki — dkd m gu — didigkm) + g np (r^rf m — T^T^) , (7) 

where the matrices T kl are the Christoffel symbols: 

r L = \a im {di9mk + d k g m i - d m g k i) . (8) 

Metrics which give zero Riemann curvature over the entire manifold are said to be flat. A 
coordinate transformation can always be performed so that any flat metric will transform 
into the Euclidean metric over the entire manifold. 

We can think of the system as a collection of particles and/or fields embedded in a 
manifold M diffeomorphic to R 3 . We will refer to M as the background manifold. M is 
a flat space and therefore the Euclidean metric can be used to compute distances over the 
entire manifold. Obviously any physically possible rearrangement is constrained to keep the 
particles embedded in M (assuming no gravitational field is present). Therefore the system 
should only be able to realize strain configurations which correspond (via the relation in Eq. 
f|) to metric tensor fields which do not produce any Riemann curvature, since M is flat. This 
implies that the virtual work relation in Eq. [5] should be rewritten to reflect this constraint. 
All flat metrics can be expressed as 

9ij = $ij + 9jUi + diUj, (9) 
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where Ui is a vector field continuous over the entire manifold. The constrained virtual work 
relation is then 



Since we are now constrained to a flat space we can choose Cartesian coordinates so that 



Notice that the divergence of the stress field has appeared in Eq. [12], so this approach 
cannot be used to uniquely specify the stress field. Therefore introducing a constraint on 
variations in the metric, has produced an ambiguity in the stress field. This is the origin of 
the aforementioned first criterion for uniqueness of the stress field. To have a unique stress 
field, the system must be allowed to realize the energetics of any applied strain configuration 
which can be represented by a continuous symmetric tensor field. 

If we hypothesize that stress fields are unique, then we must explain physically how the 
system can be altered so that the Riemann curvature at some point is non-zero. Even a 
first order change away from the Euclidean metric will generate a first order change in the 
curvature. One way which such distortions could be realized is by introducing a gravitational 
"point-like" source. Then, according to General Relativity, the system would be altered by 
the fact that the source has generated curvature within the background manifold. It is 
through this thought-experiment that the link between the definition of the stress field used 
in this work and the energy-momentum tensor in the gravitational literature is established. 
Unfortunately this is not a phenomenon that can be observed by experiment, especially on 
the quantum scale. 

In addition to the gravitational interpretation, we can also consider the meaning of non- 
Euclidean metrics within the context of the continuum theory of solids. In this theory, 
the solid within the limit of vanishing lattice parameter can be described by a vector field 




(10) 



y/9 = l: 




(11) 



(12) 
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Ui which indicates the displacement of atomic positions away from some reference (usually 
taken to be a perfect crystal)!!. If the specified atomic configuration contains defects such 
as disclinations and dislocations, then the displacement field «j will not be perfectly smooth 
and therefore will violate the integrability conditions: 

(didj - djdi) u k = 0. (13) 

Within linear elasticity theory, the information regarding the disclinations and dislocations 
be represented by a symmetric tensor field known as the defect density^: 



can 

Vij = efef n d k d m u ln , (14) 

where e^i is the Levi-Civita symbol, and Uij is the continuum strain tensor field. Notice 
that rjij is the double curl of the strain field. If rjij = over the entire manifold, then the 
strain field may be written in terms of the gradient of a displacement vector field Ui which 
obeys the condition in Eq. [13]. Specifically, 

1 



u 



'.J 



- (diUj + djUi) , (15) 
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This is known as compatibility. If there are no defects present, then the strain field is 
compatible with a smooth vector displacement field. 

The defect density has a rigorous geometric meaning. It can be shown thatS 

Vij — Rij ~ ^9ijR — Gij. (16) 

Gij is known as the Einstein curvature tensor and is related to the Riemann curvature 
through the Ricci curvature Rij = g lm Ru m j and the scalar curvature R = g^Rij. Therefore 
if a particular strain configuration in a continuum model generates Einstein curvature, it 
should be interpreted as a mathematical representation of the discontinuity in the displace- 
ment of the particles in the system due to defects. The curvature should not be thought of 
as a physical manipulation of the geometry of the background manifold. However, on the 
quantum scale, the relationship between the curvature and the compatibility of the strain 
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field requires careful consideration. It is not obvious how to assign physical meaning to, for 
example, a "defect" that is point-like in an electron gas. 

Since we must consider the energetic response of the system to non-Euclidean metrics, 
the quantum-mechanical energy must be generalized so that it is unambiguously defined 
for metrics that generate an arbitrary Riemann curvature. This is the origin of the second 
criterion for uniqueness of the quantum stress field. If we take the ab-initio perspective, 
we can build the generalized Hamiltonian from only a few facts. First it is reasonable to 
require that the energy transform as a covariant scalar regardless of the geometry. Sec- 
ond, the generalized Hamiltonian must reduce to the known correct flat-space result when 
evaluated with a metric which generates zero Riemann curvature. Finally, the Hamiltonian 
must be well-behaved in the sense that it corresponds to a quantum field theory that is 
renormalizable. Unfortunately, these conditions are not strict enough to uniquely specify a 
single Hamiltonian. There are an infinite number of possible choices, all which satisfy the 
above three conditions. In particular, the ambiguity manifests itself in the kinetic energy. 
In Euclidean space, the Hamiltonian for a spin- zero free field ip is 

H = i J V^ f • V^d 3 x. (17) 

To generalize for a manifold of arbitrary curvature, we can write H in a covariant form: 

H' = \j V99 ij di^d^d 3 x. (18) 

The operator H' satisfies the aforementioned requirements for a generalized Hamiltonian. 
Notice H' could have been constructed by writing H in terms of curvilinear coordinates in a 
flat space. This is known as the principle of least coupling. It is the simplest way to generate 
operators which are suitably generalized for curved spaces. Also it is the mathematical 
formulation of the equivalence principle in General Relativity. However, we could postulate 
another generalized form of H: 

H" = J y/g (^g ij diftdji) + aR + pR^^j + -fR ij diftdji?j d 3 x, (19) 



where a, /3, 7 are coupling constants to the various curvature terms. Note that H" is a 
scalar and gives the correct Euclidean space Hamiltonian when = 5ij. However the term 
multiplied by 7 will cause the theory to be non-renormalizable in a curved space. This 
can be determined via dimensional analysis. In the units h — 1, c — 1 (commonly used 
in relativistic field theories^) the Hamiltonian must have units of mass, or equivalently 
(length) -1 . This implies the integrand, or Hamiltonian density, must have units of (mass) 4 
since the measure d 3 x has units of (mass) -3 . (The metric is unitless.) The spin-zero field 
ip has units of (mass) 1 , which can be deduced from the first term of H" which contains 
only gradients. The scalar curvature R has units of (mass) 2 , since it contains the second 
derivatives and products of the first derivatives of the metric. Hence we can conclude that a 
has units of (mass) 2 , (3 is dimensionless, and 7 has units of (mass) -2 . If a coupling constant 
has negative mass dimension, then the theory is non-renormalizableil. Therefore we can set 
7 equal to zero. 

The ambiguity can be further reduced if we physically reason that if no matter is present 
in the system, then the stress field should be zero everywhere i.e. the part of H" multiplied 
by a is a background term. Furthermore, this term's contribution to the stress field will 
vanish when the stress field is evaluated at the Euclidean metric. Consequently we may set 
a = 0. However, we currently cannot determine any restrictions on the value of /3 from first 
principles. This is the origin of the ambiguity in the dependence of the energy on the metric. 
Therefore the quantum stress field will only be unique up to the determination of j3. We 
further note that this problem is well-known in the gravitational and constrained-dynamics 
literature, and to date remains unresolvec&~H. 

IV. THE DFT STRESS FIELD 

We now use the geometric approach to derive the quantum stress field of a many-electron 
system in the presence of a fixed set of classical positive charged ions using local density 
functional theoryilll. The ground state electronic charge density of the system is written 
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as n(x) = 0^( x )0m( x )' where M are single-particle orthonormal wavef unctions. For 
this derivation, we assume orbitals with fixed integer occupation numbers. The extension 
to metals with Fermi fillings is straightforward, simply necessitating use of the Mermin 
functional instead of the total energyS. The total charge density of the system can be 
written as a sum over all ionic charges and n: 

p(x) = £^(x-R,)-n(x), (20) 

where is the charge of the p-th ion located at position R^, and the presence of ^fg insures 
proper normalization of the delta function. The energy of the system can be written as the 
following constrained functional which is appropriately generalized for arbitrary Riemannian 
spaces: 

E = E k + ^coulomb + E xc - K (/ Vg<P*^^d 3 x - lj . (21) 

Here E k is the single particle kinetic energy including the aformentioned ambiguity, -Ecouiomb 
is the classical Coulomb interaction between the total charge density and itself, and E xc is 
the exchange-correlation energy of the electrons. The appearance of the last term in Eq. ^l| 
is due to the orthonormality constraint of the orbitals. (We choose a unitary transformation 
on {4>fj,} which enforces orthogonality.) One can express E as an integral over an energy 



densityc3E3. For convenience, we express the energy terms in Eq. 21 as the following: 

-^Coulomb = J y/g (pV - ft^^FiFj) 



j3 X 



£ xc = J Vg ne hVA (n)d s x, (22) 

where T\ = —d^V is the electric field due to the Coulomb potential V generated by p, £lda("0 
is the LDA exchange- correlation energy density, and R is the scalar curvature. To obtain 
the electronic ground-state energy, we require SE/Scj)*^ = with the additional constraints 
of a fixed metric (Sgij = 0) and a fixed ionic charge density (Sp = —Sn). This implies that 
the orbitals must obey the Euler-Lagrange equations 
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- (vqiFw*) + + ^ c ;; lomb 0, + ^„ = a^„, (23) 

which can be considered the Kohn-Sham equations for a general Riemannian manifold. Also, 
a least-action principle for -^coulomb requires that p and V obey the Poisson equation: 

(y/g^d s V) = -4vrp. (24) 

We now vary the total energy with respect to the metric. It can be proven that we do not 
need to consider variations in the electronic wavef unctions, charge density and potentials, 
since all such variations would vanish due to Eq. ^ and Eq. [23. This is the same principle 
used in the derivation of the Hellmann-Feynman force theorem and the energy-momentum 
tensor (the variation of the action with respect to metric) in general relativity@'0fi Per- 
forming the variation of the total energy with respect to the metric gives the stress field in 
local density functional theory as: 



(Tij = ~ di&ffi&v + A ij + + 9ij (7, Yl d k^d\^ 

- + ne LDA (n) - £ ^ [A„ + V^J , (25) 

where we have used the relation dy/g/dg^ = — ~y/ggij. For the moment, we simply denote 
the contribution of the kinetic ambiguity to the stress field as Aij. Using Eq. we can 
rewrite Eq. |25] as 



°a = - E di&ffi&v + A a + -^Fi^i + 9a Q E d k4>l dk 

+ A= £ W« (V99 Kk d k ^) - PRn - ^T k T k + n (e LDA (n) 



6E* 
5n 



(26) 



Eq. |26| is most useful when evaluated at the Euclidean metric. The are then solutions 
to the standard Kohn-Sham equations. From here on, we will refer to with an implied 
evaluation at the Euclidean metric. 

At this point we briefly describe how to obtain the explicit form of A^ in Euclidean 
space. We need to evaluate 
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23 r ^ (\/Vx)i?(x) . , 

My) = I Mx) w ^, ( J y) l (27) 



(Variations in n are zero due to the aforementioned use of the variational principle.) Keeping 
in mind we are only interested in the result evaluated for g^- = 6^, and Ry = 0, the only 
non-zero term is 

A rs (y) = %= [ J^n{x)g ik {x) 5 -^^d s x. (28) 



After substantial manipulation and an integration by parts, we obtain 

A rs (y) = -%L= [ J^dm^c) ( ^ fc (x)^44 - /(x)?^t1 d 3 x, (29) 



plus additional terms which are zero when g^- = 5ij. Explicit computation of the functional 
derivatives gives the final result for Euclidean space: 

A %3 = 3 (did jn - 5 l0 d k d k n) . (30) 

Therefore, the complete stress field evaluated at = is: 

ay = di&ffi&n + P (didjn - 5 ij d k d k n) + —T i J : i + ( - 0*$$^ 

+ \Y. Pn&dkh - ^ k ^k + n (e tDA (n) - ^ j . (31) 

It is important to note several key features of the form of cr^. First the ambiguity in the 
kinetic stress is identical to that identified by Godfrey via a different formalism^!. Also, the 
kinetic stress contains diagonal terms which are similar to the symmetric and antisymmetric 
kinetic energy densities. The contribution of the exchange-correlation energy to the stress 
field is only in the diagonal (pressure-like) terms, which is the proper behavior for local 
density functionals0 and is identical to the exchange-correlation stress derived in Ref. 21. 

The Coulombic contribution to cr^ is equivalent to the classical Maxwell stress field. This 
Coulombic term can be obtained by Filippetti and Fiorentini's formalism if one chooses the 
so-called Maxwell gauge0. We emphasize that it has been mathematically proven in the 
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gravitational literature that the usual energy-momentum tensor for electromagnetism is 
uniquelUll. Since this energy-momentum tensor is equivalent to the Maxwell stress in the 
non-relativistic limit and the absence of magnetic fields, the classical Coulombic interaction 
in the stress field is rigorously free of ambiguities. 

By integrating the stress field over all space, note that we can obtain the total stress T^: 

T v = J j- E diPtfirf* + ^iFj - Si—^r 1 

+ Sijn U DA (n) - j d 3 x, (32) 

which is identical to the expression derived by Nielsen and Martin0. 

It is also instructive to consider the form of the pressure field P defined as one-third the 
trace of the stress field. Using our formulation for we obtain: 



P = I E P^Mp ~ ^ k d k n - + \ E <Pld k dk<P» + n ^ LDA ( 



-) - S ~t) • (33) 



Notice the pressure field contains the Maxwell energy density and a combination of the sym- 
metric and antisymmetric kinetic energy densities. The contribution of the LDA exchange- 
correlation pressure is consistent with the unambiguous exchange-correlation energy density 
derived by Burke and co-workers0. 

Also present is the undetermined contribution of the Laplacian of the charge density. We 
should consider whether this kinetic ambiguity impairs us from using the the pressure field 
to obtain novel qualitative information from total energy calculations. We would therefore 
need to determine a realistic range for the parameter f3. Such a range can be derived if 
we consider the various forms of the kinetic energy density. Symmetric and antisymmetric 
forms can be converted from one to the other by adding a term proportional to the Laplacian 
of the charge density: 



~ E d k <p;d k ^ - -d k d k n = -W. + <P,d k d k <p;) . (34) 

Hence a realistic range for (3 in Eq. |33] would be from to 1/8. 
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V. APPLICATION: SOLID HYDROGEN 



In order to demonstrate the utility of our stress field formalism, we compute within DFT 
the pressure field, for two phases of solid molecular hydrogen under external hydrostatic 
pressure of 50 GPa.EH. Both structures consist of stacked two-dimensional triangular lattices 
of hydrogen molecules, with the molecular axis parallel to the stacking direction and a repeat 
unit of two layers. The m-hcp structure has alternating layers shifted so that each hydrogen 
molecule is directly above triangular hollow sites in the neighboring layers. The second 
structure (belonging to the Cmca space group) has a different shift, so that each molecule lies 
directly above midpoints between nearest-neighbor pairs of molecules in adjacent layers@. 
The energetics and electronic properties of both structures have been studied extensively 
from first principles^llii. 

Examination of the pressure field permits us to rationalize the energy ordering of these 
structures. The m-hcp structure is energetically favored by 60 meV/molecule. Figure |]A 
shows a contour plot of the pressure field for the Cmca structure. Here we have set the 
constant j3 equal to zero. The pressure is tensile (greater than zero) through the interstitial 
region, indicating that contraction is locally favorable. The tensile pressure is greatest in 
the volume directly above and below each molecule, averaging 3 eV/A 3 . This implies that 
the system would energetically favor increased intermolecular coordination. In Figure [TJB 
we show a similar plot for the m-hcp structure. Again the pressure within the interstitial 
region is tensile. However, the pressure field has significantly rearranged, and the pressures 
in the regions above and below molecules have been reduced to approximately 2.25 eV/A 3 . 
It is also clear from the charge density plots (Figure |lp and |lp) that the reduction in 
pressure is correlated with an increase in bonding between molecules and increased charge 
derealization. The pressure fields within the molecules are compressive, and they are several 
orders of magnitude larger than the interstitial features. However, because these large fields 
are very similar in both structural phases (and the free molecule), they are not important for 
understanding relative phase stability. Thus, changes in the pressure stress field highlight 
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regions and charge density features that contribute to favorable energetic changes. 

Figure ^ contrasts the pressure fields for (3 = and (3 = 1/8. Note that there is no 
qualitative change in the pressure field within the interstitial regions since the charge density 
in this area is quite smooth and hence the value of its Laplacian is small (averaging 0.05 
eV/A 3 ). Therefore our above analysis holds regardless of our choice of (3. 

VI. CONCLUSION 

We have developed a new formulation of the quantum stress field using differential ge- 
ometry. We identify the stress tensor field as a variation of the total ground state quantum- 
mechanical energy with respect to the metric tensor field. Within this formulation we have 
uniquely specified the stress field up to a single ambiguity. The origin of the ambiguity 
lies in the inability to uniquely describe the dependence of the kinetic energy on the metric 
tensor. 

We have successfully extended this formulation to local density functional theory. Terms 
identical to the Maxwell stress can be found within the LDA-DFT stress field. We note that 
within the gravitational literature, it is well-known that the Maxwell stress is unique. We 
describe the similarities between the pressure field obtained by our method and the various 
terms found in the total energy density. We then show that the pressure field contains an 
ambiguity in the form of the Laplacian of the charge density multiplied by a constant (3. 
Through the relationship between terms in the pressure field and the kinetic energy density 
we offer a realistic range of choices for the value of (3. 

To demonstrate the usefulness of stress field as a qualitative tool in combination with 
total energy calculations, we used the pressure field to rationalize the energetics of two 
phases of solid molecular hydrogen. We finally show that our conclusions remain invariant 
with respect to realistic choices of (3. 
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FIGURES 



FIG. 1. Contour plots of the pressure field with (3 = and the charge density within DFT for 
the Cmca structure (panels A and C) and for the m-hcp structure (B and D) of solid hydrogen. 
The vertical axis is the stacking direction for the layers, and the horizontal axis is the direction 
along which alternating layers are shifted. The plots are 6.794 A high and 7.206 A wide. Ten 
contours are shown, over a range of 0-3 eV/A 3 for the pressure fields, and from 0-2.5 e/A 3 for the 
charge densities. 

FIG. 2. Contour plots of the pressure field in solid hydrogen for the Cmca and m-hcp structures 
with (3 = (panels A and C), and (3 = 1/8 (panels B and D). The vertical axis is the stacking 
direction for the layers, and the horizontal axis is the direction along which alternating layers are 
shifted. The plots are 6.794 A high and 7.206 A wide. Ten contours are shown, over a range of 
0-3 eV/A 3 . 
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